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Abstract: The drought has enormous adverse effects on agriculture, water resources and environment, and 
causes damages around the world. Drought risk assessment and prioritization of drought management can 
help decision makers and planners to manage the adverse effects of drought. This paper aims to determine 
the risk of drought in Iran. At the first stage, standardized precipitation index (SPI) was calculated for the 
period 1981—2016. Then the probability map of different drought classes or drought hazard probability 
map were prepared. After that the indicator-based vulnerability assessment method was used to determine 
the drought vulnerability index. Five indices including climate, topography, waterway density, land use and 
groundwater resources were chosen as the most critical factors of drought in Iran and followed by the 
analytical hierarchy process questionnaire, the weights of each index were obtained based on expert 
opinions. Fuzzy membership maps of each index and sub-index were prepared using ArcGIS software. The 
drought vulnerability map of Iran was plotted using these weights and maps of each indicator. Finally, the 
drought risk map of Iran was provided by multiplying drought hazard and vulnerability maps. According to 
the 43-completed questionnaires by experts, climate index has the highest vulnerability to drought. Climate 
does not have an important role in drought hazard index, but it is the most crucial factor to classified 
drought vulnerability index. The results showed that central, northeast, southeast and west parts of Iran are 
at high risks of drought. There are regions with different risks in Iran due to unusual weather and climatic 
conditions. We realized that the climate and the groundwater situation is almost the same in the central, east 
and south parts of Iran, because the land use plays a crucial role in the drought vulnerability and risk in 
these areas. The drought risk decreases from the center of Iran to the southwest and northwest. 
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1 Introduction 


Drought is one of the most destructive natural phenomena which have devastating impacts on 
agriculture, water resources and environment (Dow, 2010; Popova et al., 2014; Yu et al., 2014). 
This phenomenon starts slowly but can take a long time and covers a large area compared to the 
other natural hazards (Rossi et al., 1992; Mishra and Singh, 2010; Lena et al., 2014). In the future, 
with an increase in severe weather incidents, it is expected that drought will occur more often and 
stronger (Trenberth et al., 2014). 

The definitions of drought are classified as meteorological drought, hydrologic drought, 
agronomic and social economic drought (Dracup et al., 1980; Orville, 1990; Wu et al., 2016). 
Meteorological drought is the most important type (Nasrollahi et al., 2018), which is the driving 
force behind other droughts, resulting in a reduction in the average long-term precipitation. With 
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the meteorological drought, its adverse effects appear and the other types of drought will gradually 
take place and the combind adverse impacts can cause the destructive impacts in an area. 

Economically speaking, drought is the most dangerous natural phenomenon during history 
(Karabulut, 2015) and has always endangered global food security (Dai, 2011; Mu et al., 2013). 
Owing to the consumption of water resources in recent decades, new aspects of drought has 
emerged, and it has affected the enormous human population around the world for years (Kim et 
al., 2013; Lei et al., 2016). Hence, how to manage drought is one of the essential challenges 
worldwide (He et al., 2013). 

In recent years, a risk-oriented approach has been adopted for disaster risk management (Tonini 
et al., 2012; Skakun et al., 2014). In this regard, governments have attempted to develop a drought 
management policy with emphasizing on risk management instead of the traditional approach 
(Wilhite, 2000; Sivakumar and Wilhite, 2002). Evaluation of drought risk is a two-argument 
function including drought hazard and vulnerability (Jonkman et al., 2003; IPCC, 2007; Jayanthi 
et al., 2013). 

Drought hazard assessment is the vital aspect of drought risk (Hao et al., 2012), which means 
investigation of probability for drought intensity and its frequency (Zhang and Zhang, 2016). 
Repetition of high-intensity drought can have hazardous effects (IPCC, 2007). Drought hazard 
describes physical aspects of drought and has an essential role in the connection between drought 
risk and vulnerability (He et al., 2011). Therefore, drought hazard is the most important aspect of 
drought risk (Hao et al., 2012). Drought hazard can be defined as a probability for drought intensity 
and its frequency (Zhang and Zhang, 2016). In other words, drought hazard is the probability of 
drought occurrence. The drought hazard index is calculated by statistical analysis of drought indices 
like standardized precipitation index (SPI), standardized precipitation evapotranspiration index 
(SPEI) and Palmer drought severity index (Zargar et al., 2011; Beven et al., 2018). 

The vulnerability is a broad concept with various definitions (Easter, 1999; Birkmann, 2006; 
Fiissel, 2007). Vulnerability means human and natural potential sensitivity of an area when drought 
occurs. In other words, the differences between natural properties and economical-social features 
caused different sensitivity when drought occurs (Santos et al., 2001). There are different ways to 
measure drought vulnerability. One of the most common methods is index-based vulnerability 
assessment (IBVA) methods. To quantify the vulnerability, IBVA depends on several potential 
indices and combining them with the related levels of vulnerability. Many researchers have used 
this method in local and global scales (Moss et al., 2001; Leichenko and O-Brien, 2002; Vincent, 
2004; Ashok and Sasikala, 2012). There are two options for calculating the vulnerability using 
IBVA. The first is to consider the equal importance and therefore equal weight for all factors (Me- 
Bar and Valdez, 2005); and the second is to find the different emphasis and to assign different 
weights to the factors causing vulnerability. In the second option, there are different approaches to 
consider different weights. To determine the weights of the factors influencing vulnerability, 
techniques such as the analytical hierarchy process (AHP) can be used. 

Many studies on drought hazard, vulnerability to drought and drought risk have been carried out. 
For example, Sahahid and Behravan (2008) used several physical factors such as soil available 
water capacity, production amount and percentage of the irrigated area; and socio-economic factors 
such as population density, gender ration and the rate of people employed in the agriculture sector 
to provide the map of vulnerability to drought in the west of Bangladesh. 

Han et al. (2016) calculated the comprehensive drought risk index of southwestern China using 
climatic, geographic, soil and remote sensing data, with the help of AHP. The results of the research 
indicated that prevention capacity factor and reducing drought, drought hazard, sensitivity and 
environmental vulnerability are involved in calculate the drought risk index. 

Nasrollahi et al. (2018) studied the drought risk of Semnan province of Iran. In this research 
drought hazard was calculated using SPI. Vulnerability of the area to drought was evaluated by 
indices such as agricultural land and water used in domestic, agronomy and industry sections. The 
results show that the most significant part of the study area has the medium (62.46%) and very 
intense (53.12%) drought risk. Chen et al. (2019) evaluated drought risk in agriculture in China 
using SPI to determine the drought hazard. Results show that agricultural drought risk decrease 
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from northeast to southwest. They suggested that the research results can be used to mitigate the 
effects of future droughts in China. To study the future drought risk in Africa, Ahmadalipour et al. 
(2019) calculated drought hazard index using SPEI. The drought vulnerability index was calculated 
by available layer such as land use and water resources, etc. The drought risk was predicted based 
on two climate scenarios for 2100. The results indicate that Africa would face an unprecedented 
drought hazard if there is no risk management plan against to climate change. 

Drought is one of the most costly natural disasters (Lee et al., 1999; Khoshnodifar et al., 2012), 
having adverse environmental, social and economic impacts in Iran. The researches on climate 
change indicate that the frequency and duration of drought in the Middle East will also increase 
(Dai, 2011). Government activities in Iran have focused on drought crisis management (response 
to drought and recovery), but supervision and increase readiness to copy with drought as risk 
management has not been considered (Nasrollahi et al., 2018). So, assessment of drought risk and 
prioritizing management strategies to cope with this phenomenon is necessary. Therefore this paper 
tries to demonstrate the probability of drought map of Iran to determine drought-vulnerability areas 
in Iran for management when severe drought occurs. 


2 Materials and methods 


2.1 Study area 


Iran (1.65x10° km’; 25°03'-39°47'N, 44°14'-63°20’E) has very different weather conditions due 
to its broad ecological diversity and terrain variety (Heydari et al., 2018; Behrang, 2019). Figure 1 
shows the climatic diversity of Iran, which was prepared by mixing different classes of climate map 
developed by Iran Forests, Range and Watershed Management Organization 
(http://www. frw.org.ir/00/En/default.aspx). 
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Fig. 1 Location of the study region, distribution of synoptic stations and climatic diversity in Iran 


2.2 Data 

The data of 46 meteorological stations including the mean annual precipitation, temperature and 
evaporation were obtained from the Iran Meteorological Organization (http:// http://www.irimo.ir). 
The meteorological data were interpolated using the IDW (inverse distance weighting) method in 
ArcGIS 10.3 (Esri, USA). Locations of the meteorological and interpolated stations in Iran are 
shown in Figure 1. For groundwater resources condition, the data of piezometer wells prepared by 
Iran Water Resources Management Company (https://www.wrm.it/) were used. The waterways 
layer of Iran was prepared by Forests Range and Watershed Management Organization of Iran 
(http://www.frw.org.ir/02/En/). The land use layer of Iran was provided by Forests Range and 
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Watershed Management Organization of Iran (http://www.frw.org.it/02/En/). For topography index 
the DEM (digitial elevation model) 90 m of Iran (https://cgiarcsi.community/data/srtm-90m- 
digital-elevation-database-v4-1/) was used to provide slope, aspect and height layers. 


3 Methodologies 


3.1 Drought hazard 


The monthly precipitation data of 46 synoptic stations with long-term (1981—2016) and reliable 
statistics were selected to calculate the drought hazard (Fig. 1). The data of station were completed 
and their missing data were estimated from the nearest neighbor stations using linear regression. 
Then, the SPI was used to determine the drought hazard (Potop et al., 2012). 

In order to calculate the SPI, the precipitation data of the statistical period of each station were 
fitted to the gamma probability distribution function, g(x), as Equation 1. 
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where a is the shape parameter, which is estimated using the method of maximum likelihood; £ is 


the scale parameter; x is the precipitation (mm); and T (æ) is characteristics of the gamma function. 
The parameters of g(x) were calculated by the following equations: 
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where A is a dimensionless component used for calculation of a and J; n is the number of rain 
observations; and x is the average precipitation (mm) for a given month during the statistical 
period. The calculated parameters allow the distribution of precipitation at the station be effectively 
represented by a cumulative mathematical probability function, G(x), as Equation 5. 


x 1 x i 
G@) = f - f xt eB dx. (5) 
3 BV (a), 
Since the gamma function is not defined for x=0, the distribution of the precipitation become 
zero. 


A(x)=q+U-q)G(x), (6) 


where H(x) is the cumulative probability function if x=0 and q is zero probability. H(x) then 
converted to the standard normal distribution for calculating SPI values ((Mckee et al., 1993). 

To simplify operations and increase the accuracy of calculations, the DIP (Drought Index 
Package) software version 2 was used to calculate SPI in a 12-month time step (Morid et al., 2005). 
After calculating the SPI, the classification of the drought severity was performed (Table 1). 

In this study, drought hazard was evaluated by the probability of drought occurrence. To calculate 
the drought probability, at first, the occurrence frequency of different drought classes including 
moderate, severe and too severe was determined in each synoptic station. Then, by a division of 
probability of drought occurrence in each category to all possible states, the percentage of drought 
probability in each category was calculated. In the next step, the point data interpolation was 
performed using geostatistical analysis tools through IDW method in ArcGIS 10.3. 


Table 1 SPI (standardized precipitation index) classification (Mckee et al., 1993) 
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Drought category SPI value 
Too severe drought < -2.00 
Severe drought —1.99<SPI< —1.50 
Medium drought —1.49<SPI< —1.00 
Normal —0.99<SPI<0.99 
Medium wet 1.00<SPI<1.49 
Severe wet 1.50<SPI<2.00 
Too severe wet >2.00 


Because various drought intensities do not have the same value, for quantifying and mapping the 
drought hazard, the weight was given to each class of drought according to Table 1. The rating of 
1 to 4 was assigned to the moderate, severe and very severe drought groups (Table 2). The map of 
each intensity class was classified into four subclasses based on percentage probability value using 
natural break method. This method based on a computational algorithm tries to minimize the 
differences between the data in each class and to maximize the distinctions among categories. This 
algorithm uses the average of each range for creating classes to achieve greater uniformity for 
distributing data. In this study, ArcGIS 10.3 was used to determine natural break points. Then, 
number 1 and 4 were assigned to the subclass with the lowest and highest probabilities of drought, 
respectively. 

Finally, the drought hazard index (DHI) was calculated using Equation 7 and ArcGIS 10.3 in a 
time step of 12 months: 


DHI = (MD; x MDw) + (SDr x SDw) + (VSDr x VSDw) , (7) 
where MD,, SD,, and VSD; are given rating to the moderate drought (MD), severe drought (SD) 
and very severe drought (VSD), respectively, and MDw, SDw, and VSDw are weights for moderate, 


severe and very severe drought. According to Equation 1, DHI can have a value of 6-24. Then, 
DHI maps can be classified into four classes of low, moderate, severe and very severe. 


Table2 Weight and rating assigned to the drought in a time step of 12 months 


Drought class Weight Rating Occurrence probability” (%) Area (%) 
1 <6.96 20.02 
2 6.96-8.69 42.23 
Moderate drought 1 
3 8.69-10.56 28.89 
4 >10.56 8.86 
1 <3.76 32.21 
2 3.76-5.05 39.92 
Severe drought 2 
3 5.05-6.84 14.89 
4 >6.84 12.98 
1 <1.59 12.47 
2 1.59-2.34 35.43 
Very severe drought 3 
3 2.34-3.03 43.18 
4 >3.03 8.92 


Note: *, the occurrence probability of each rating was calculated using natural break method in ArcGIS 10.3. 


3.2 Drought vulnerability 


In this study, IB VA method was used to calculate the drought vulnerability. At first, effective indices 
of drought vulnerability was determined. They are selected according to the purpose of the study 
and properties of the study areas (Hinkel, 2011; Araya-Mufioz et al., 2017) and their variation 
(Dabanli, 2018). Based on the previous studies and consultation with Iranian experts in the field of 
drought, five indices including climate, topography, waterway density, land use and groundwater 
resources were selected and for each index some sub-indices were considered (Fig. 2). 

The AHP method, e.g., pairwise comparison matrix, was used for weighing each effective index 
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and sub-indices (Sadeghravesh et al., 2015). After the formation of the pairwise matrix, if the matrix 
contrariety rate was acceptable, the weight of each criterion was calculated (Malczewski, 1999; 
Hill et al., 2005). The 43-completed questionnaires by drought experts were used for analyzing and 
estimating the weight of each index and sub-index. Also, Expert Choice (Expert Choice Inc., USA) 
was used for formation of the pairwise matrix, calculating the contrariety rate, relative and final 
weights of criteria. 


Depletion 


Drought vulnerability 
Wat 
f oe, Climate Land use 
ensity 


Fig. 2 Analytic hierarchy process (AHP) for drought vulnerability determination 


Groundwater 
resources 


Topography 


In the next step, the map of each index and sub-index were prepared in ArcGIS 10.3. The selected 
indices, sub-indices and the method for calculating them are described as follows: 


Groundwater resources index: groundwater is one of the essential factors which affects 
drought vulnerability (Jordaan et al., 2019). In this study, groundwater resources index was 
divided into two sub-indices: water depth average and annual groundwater drop level. The 
fundamental data of these two indices were calculated by the reported 17-a data (2000- 
2016). A special water depth average and an average amount of annual groundwater drop 
level were attributed to each sub-watershed. Maximum weight was for the watershed which 
had the highest amount of drop and depth; and this watershed was much sensitive to drought 
and vice versa. 

Waterway density index: high density of waterway shows much runoff water. Consequently, 
much runoff in an area increases sensitivity against drought (Jain et al., 2015). 

Climate index: climate is one of the most important factors in determining drought 
vulnerability. The index of climate included three sub-indices of precipitation, evaporation 
and temperature. For the preparation layer of these three sub-indices, data of mean 
precipitation, evaporation and annual temperature of 46 synoptic stations of the country in 
the period 1981—2016 were used. Then, zoning was conducted for the three sub-indices using 
conventional Kriging method in ArcGIS 10.3. Increasing precipitation decreases drought 
severity, so regions with more precipitation will have a lower weight value. Raising 
temperature increases the susceptibility to drought, accordingly, the maximum weight value 
assigned to the areas with the maximum temperature, and the minimum weight value was 
allocated to areas with the minimum temperature. The difference between evaporation and 
precipitation can be considered as a factor in assessing drought. Increasing this difference 
will increase the sensitivity to drought. Thus, areas with more difference allocated more 
weight. 

Land use index: one of the most important indices which affect vulnerability is land use 
(Nelson et al., 2005; Dayal et al., 2018). The layer of the index was prepared using ArcGIS 
10.3 and utilizing land use layer. At the next step, the maximum and minimum weight were 
given to the areaes with irrigated cultivation and rangeland or bare lands respectively. 
Topography index: since the topography has an enormous effect on available water, it could 
be one of important factors for drought vulnerability. This index are composed of three sub- 
indices of slope, aspect and height. As we know, increasing height reduces sensitivity to 
drought, so the areas with the maximum height assigned less value and the areas with the 
lowest one assigned more value. The increased slope leads to water being unreachable to the 
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plants rapidly; thus, areas with greater slopes were assigned higher weight value. In the 
standardized map of direction parameter, the highest weight attached to south, west, east and 
north slopes, respectively. 

Finally, the fuzzy logic was used to combine and overlay the layers of indices and sub-indices. 
In fuzzy logic, membership of an element in a set is defined with a value in the range of 0 to 1 (Sui, 
1992; Bonham-Carter, 1994). Each pixel of indices and sub-indices maps was gotten the value 
between 0 and 1, and O indicates a little impact while 1 suggests the high impact on drought 
vulnerability. 

Fuzzy operation gamma was resulted by multiplying operators of fuzzy algebraic product (Eq. 8). 


u(x) = (Fuzzy algebric sum)” x (Fuzzy algebric production)!” , (8) 


where u(x) represents the fuzzy operation gamma equation and y is a number between 0 and 1. The 
correct and conscious selection of y creates some values in outputs that shows flexible adaptation 
between downward tendencies of fuzzy multiplication and upward trends of fuzzy summation 
(Bonham-Carter, 1994). In this study, y was considered as 0.9. Dayal et al. (2018) used y=0.9 to 
assess drought risk in southeastern Australia and achieved acceptable results. The map of drought 
vulnerability index was obtained by overlaying indices and sub-indices maps. 


3.3 Drought risk 


By calculating drought hazard index (DHI) and drought vulnerability index (DVT), the drought risk 
index (DRI; Eq. 9) was calculated (Knutson et al., 1998; Downing et al., 1999; Downing and 
Bakker, 2000; Wilhite, 2000). 


DRI=DHIx DVI. (9) 


Drought risk will increase with DRI. In the end, each map (drought hazard, vulnerability and 
risk maps) was classified into five classes as very low, low, moderate, high and very high by 
utilizing natural break in ArcGIS 10.3. 


4 Results 


4.1 Drought hazard 


The 12-month severe drought percentage probability map (Fig. 3) shows that the northeastern, 
central, southwestern, and part of northwestern Iran have high drought probability. The coasts near 
the Strait of Hormuz, southern and the north coast of Iran and around the lake of Urmia have a very 
low probability of such a severe drought for 12 months. In general, the eastern and southeastern 
Iran have the highest possibility of severe drought and the northwest, northeast coast of the Persian 
Gulf to the center of Iran have the minimum probability. It can be concluded that the dispersed 
parts of northwest to southeast have the highest probability, while the west, east and northeast parts 
of Iran have the lowest drought probability. 


4.2 Drought vulnerability 


The 12-month DHI spatial distribution in Iran for the period 1981-2016 (Fig. 4) shows that the 
middle belt of Iran from east to west and the southwest have high probability of drought in a time 
step of 12 months. The coast of the Strait of Hormuz to the center of Fars province and the northwest 
region has the lowest hazard of drought for the period. Almost half of Iran located in the very high 
class (48.2%) and the smallest area are also related to the low class (9.9%). More than 65% of the 
country is in high and very high classes. 

The weights of each index and sub-index were calculated by AHP to determine the drought 
vulnerability. Table 3 shows that the climate index has the highest susceptibility to drought. In other 
words, climate index plays a major role in deciding the vulnerable areas (more than 50%). The sub- 
indices of precipitation, slope and depletion have the most significant impacts on climate, 
topography and groundwater resources, respectively. 
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Fig.3 Drought probability map for 12 months: moderate (a); severe (b); and very severe (c) 
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Fig. 4 Drought hazard index for time step of 12 months 


Table 3 Index and sub-index weights determined by the analytic hierarchy process 


202101.00053v1 


chinaXiv 


ChinaX ivé EAT! 


JOURNAL OF ARID LAND 
Index Weight Sub-index Weight 

Precipitation 0.693 
Climate 0.516 Evaporation 0.088 
Temperature 0.219 

Waterway density 0.040 - - 
Slope 0.490 
Topography 0.061 Aspect 0.197 
Height 0.313 

Land use 0.164 - - 
Depletion 0.875 

Groundwater resources 0.219 

Depth 0.125 


The fuzzy membership map of waterway density index (Fig. 5) shows that some regions in north- 
central and eastern Iran are less vulnerable to drought, and southwest and western Iran are more 
susceptible to drought due to the high density of waterway. 
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Fig. 5 Fuzzy membership map of waterway density index 
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Fuzzy membership maps of land use index (Fig. 6) shows that areas with more vegetation cover 
have a more significant role in determining vulnerability to drought. The central Iran, especially 
the Dasht-e Kavir and Kavir-e-Lut, have no critical role concerning land use and vulnerability to 
drought. The fuzzy membership map of groundwater resources index (Fig. 7) shows that each 
watershed has its requirements depending on the management of groundwater resources in that 
watershed. Fuzzy membership topographic map index does not show a specific trend, and each 
pixel has been weighed according to its slope, direction and height (Fig. 8). Checking the fuzzy 
membership map of climate index shows that areas in the central Iran are more sensitive to drought 
and are more vulnerable compared to the western, northwestern and northern regions (Fig. 9). After 
determining fuzzy membership of indices, the map of vulnerability to the drought of Iran was drawn 
using fuzzy overlap in ArcGIS 10.3, and was classified into five classes (Fig. 10). 

Based on the drought vulnerability map of Iran, the central, northeastern and southeastern Iran 
are located in very high class, but southeastern, eastern and north-central part are located in very 
low class. Highlands of Zagros and Alborz mountains are also found in high class and most parts 
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of the northwestern and western Iran as well as the north coast were in the moderate and low classes. 
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Fig. 7 Fuzzy membership map of groundwater index (a) and its sub-indices: groundwater depletion (b) and 
groundwater depth (c) 
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Fig. 8 Fuzzy membership map of topography index (a) and its sub-indices: aspect (b), height (c), and slope (d) 
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Fig. 9 Fuzzy membership map of climate index (a) and its sub-indices: precipitation (b), temperature (c) and 
evaporation (d) 
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Fig. 10 Drought vulnerability map of Iran 
4.3 Drought risk 


Drought risk map of Iran was prepared by multiplying the two sheets of hazard and vulnerability 
to drought, which was classified into five classes (Fig. 11). The results showed that parts of central, 
northeast, southeast and west parts of Iran are at high risk drought class. North and northwest parts 
of Iran are at low and very low risk classes. 
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Fig. 11 Drought risk map of Iran 
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5 Discussion 


Drought is a complicated phenomenon and comprehensive drought risk management is necessary 
to reduce its negative effects. In this study, the DRI of Iran was evaluated using DHI and DVI. DHI 
was estimated from a study of occurrence probability with different drought intensity, which is 
calculated by SPI. In other words, drought occurrence probability has a major effect on DHI (Wang 
et al., 2012; Liu et al., 2014; Han et al., 2016). In general, the results show that, by increasing 
drought intensity, the probability of drought occurrence will decrease. The research results of Li et 
al. (2018) also proved this. Drought probability maps (Fig. 3) indicate that drought occurrence with 
different intensities is not a function of climate in different regions. Therefore, DHI is not dependent 
on climate. It is just because drought will occur in different types of climate (Vicente-Serrano et 
al., 2012; Wang et al., 2017; Lai et al., 2019) and the intensity of the drought will not relate to 
climate. Also, Iran's DHI indicates that an area located in very high risk class might have lower 
probability of very severe drought risks than an area located in a low risk class. For example, the 
probability of 12-month very severe drought is low in eastern Iran, but this area is located in the 
very high drought risk class. This indicates that DHI is a combination of different drought 
intensities and can be used for drought-related planning. 

The assessment of drought risk is not possible without determining drought vulnerability (Jain 
et al., 2015; Ekrami et al., 2016; Thomas et al., 2016). In this study, IBVA method was used to 
determine drought vulnerability. Five indices and their sub-indices were considered as effective 
factors on drought vulnerability. There are a lot of methods to determine the weight of indices, but 
it is notable that determining the precise weight of indices is not easy (Kiem et al., 2016; Tanago 
et al., 2016; Muyambo et al., 2017). We used AHP method to determine the weights. Table 3 shows 
that the climate index has the most significant weight and has the most efficacy on drought 
vulnerability. Figure 9 also indicates that, if the climate of an area is dry, it would be more 
vulnerable to drought. 

The groundwater resources and land use are also important indices effecting DVI and have direct 
relations with human activities, which can reflect the economic and social conditions in Iran. The 
main source of water supply in Iran is groundwater (Jafary and Bradley, 2018). Groundwater 
resources has a vital role in drought vulnerability because of the high dependency in the agriculture 
sector in Iran. Non-systematic and uncontrolled consumption caused to face a crisis of groundwater 
(Hojjati and Boustani, 2010; Izady et al., 2012). Hence, crisis and risk management of groundwater 
in Iran has a vital role in vulnerability and drought risk. The significant influence of land use on 
determining DVI demonstrated that time, period and drought intensity have different effects on 
different land uses (Khosravi et al., 2017). 

The map of drought vulnerability (Fig. 10) shows that the regions nearby Urmia Lake located in 
moderate and high vulnerable classes have to be considered for drought risk management. The 
effects of drought and human activities, especially agriculture and utilizing groundwater, have 
caused a decrease in the lake water level (Ghale et al., 2018). West coasts of the Caspian Sea have 
less vulnerability against drought in comparison with the east shores because of the humid climate 
of the west. Based on selected indices, Alborz and Zagros Mountains are located in high drought 
vulnerable class. 

The evaluation results of drought risk in Iran show that risks vary in different areas owing to the 
different climatic conditions. Most of the areas in semi-arid to humid climate (Fig. 10) are classified 
as high to very high drought risk. Most of the Iranian population are situated in these areas. Hence, 
the drought risk management in these areas has special importance. 

Having hyper-arid and arid climate (Fig. 1) as well as the same groundwater condition, land use 
index has the decisive role on drought vulnerability and risk in the central, eastern and southern 
Iran. For instance, most arid areas, like Lut and central desert with high and very high drought 
probability are located in low drought risk class. From the central part of Iran to the southwest, 
drought risk gradually decreases, whereas the very high class of drought risk increases from 
northwest to Iran's central region. Land use is a sign of interaction between human and nature 
(Moghadasi et al., 2018). Therefore, human activities excessive exploitation of natural resources 
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have a decisive role in management drought risk in dry and very dry parts of Iran. 


6 Conclusions 


In this research, we tried to provide the Iran's drought hazard, vulnerability and risk maps to 
determine drought vulnerability areas when severe and very severe drought occurs. The monthly 
precipitation of 46 synoptic stations with 35-a data were used to calculate the SPI in a 12-month 
time step. Drought hazard map was provided by the probability map of drought occurrence. Based 
on the IBVA method, five indices including groundwater resources, waterway density, climate, land 
use and topography were selected as the most drought-sensitive indices in Iran to prepare the 
drought vulnerability map. Multiplying drought hazard and vulnerability maps, the drought risk 
map of Iran was created. The drought hazard map in12-month time step showed that more than 
60% of Iran have high and very high probability of drought hazard. The eastern and southeastern 
Iran have the highest possibility of severe drought and the northwest, northeast coast of the Persian 
Gulf to the central part of Iran have the minimum probability. Based on the drought vulnerability 
map, the central, northeast and southeast parts of Iran located in very high class are the most 
sensitive area to drought. This is due to the topography and waterway density. Among the sub- 
indices of topography, the weight of slope was greater than the height and aspect. In arid and hyper- 
arid regions of Iran, slope changes are much less than those in the semi-arid and humid areas, so 
the weight of the topography was low in the arid and hyper-arid regions. The waterway density is 
low in the arid and arid regions, so this index also has a lower weight than that in the semi-arid and 
humid areas. 

The results of this study illustrate that climate, groundwater resources and land use have 
important impacts on drought vulnerability in Iran. The groundwater resources and land use as 
human activities show the human role on drought vulnerability and risk. Due to population density 
and its pressure on resources, the weights of land use and groundwater resources in semi-arid and 
humid areas were greater than those in the arid and hyper-arid areas. 

Based on the drought hazard and vulnerability maps, the climate index does not have an essential 
role on drought hazard index, but it is the most crucial factor for determining the drought 
vulnerability index. The drought risk map showed that the central part, northeast, southeast and 
west areas of Iran are at high risk drought class. Some reducing-drought risk activities, including 
providing a drought monitoring and foresight system, developing the jobs not depending on the 
water, and training society to create a new method, are proposed to improve water consumption 
efficiency in industrial, agriculture, and domestic sectors. 

Drought risk map of Iran illustrated that in the dry and very dry area, land use has a unique role 
on drought risk. In other words, in these area, there is a close relationship between land use and 
drought risk. Increasing drought risk can change land use; and the land use management can effect 
drought risk and vulnerability. The hyper-arid and arid regions having a high probability of drought 
are located in very low and low drought risk classes. These areas have less management priority. 
In general, natural risk management is more important in places where human interests are at stake. 
In arid and hyper-arid regions, the population density is lower than that in the semi-arid and humid 
regions, so human resources are less pressured and land use is more stable. The combination of 
these reasons led to a lower risk of drought in arid and arid regions than in humid and semi-arid 
regions. 

It is recommended that the risk of drought in each land use be separately assessed under different 
topographic climatic conditions. This provides better risk management in arid and arid regions, and 
reduces the importance of population centers in determining vulnerable areas. Furthermore, it is 
suggested that vegetation changes as a new index in assessing drought vulnerability be achieved 
using satellite imagery. 
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